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ABSTRACT 

The Blandford-Znajek process, one of the most promising model for powering the 
relativistic jets from black holes, was initially introduced as a mechanism in which the 
magnetic fields extract energy from a rotating black hole. We study the evolution of 
force-free electromagnetic fields on regular spacetimes with an ergosphere, which are 
generated by rapidly rotating stars. Our conclusive results confirm previous works, 
claiming that the Blandford-Znajek mechanism is not directly related to the horizon 
of the black hole. We also show that the radiated energy depends exponentially on the 
compactness of the star. 
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1 INTRODUCTION 

The Blandford-Znajek (BZ) process is one of the leading 
models to explain the launching of powerful relativistic jets 
emerging from the supermassive black holes at the center 
of the galaxies (i.e. Active Galactic Nuclei), and the more 
moderated ones coming from stellar mass black holes (i.e. 
microquasars) . The main ingredients of this process are a 
central rotating black hole and an accretion disk, which sup- 
ports a magnetic field threading the black hole horizon. This 
magnetic field is twisted by the spinning black hole, produc- 
ing an outgoing electromagnetic flux which extracts energy 
and angular momentum from the spacetime. 

Although the BZ model was introduce d a long time ago 
(|R. D. Blandford and R. L. Znaiekl (|l977l 'l'l. it is only re- 
cently that many issues and theoretical discoveries concern- 
ing this mechanism have been settled. These advances on 
the understanding of the BZ process have been enabled by 
numerical simulations. For instance, it has been shown that 
only the magnetic fields lines threading the ergosphere of the 
black hole (i.e. the region near the black hole where negative 
killing energies can exist) rotate due to the fra me dragging 
effect, whether or not t hey cross the horizon (|Komissarovl 
l|2002l . |2004 l2005l . 120091 ') ) . These twisted magnetic fields are 
carrying the energy of the relativistic jet, which seems to 
come from the ergosphere. Moreover, it is now understood 
the dependence of the luminosity on the black hole spin 
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magni tude (jTchekhovskov et al.l (120101): iPalenzuela et al.l 



|2010| )) and its orientation jPalenzuela et al.l (|2010l ')'). It 



has also been shown the robustness of the process with 
respect to different boundary conditions (jPalenzuela et al.l 
(2011)), and its resemblance to ideal MHD solution s in th e 
limit of high magn etization (JMcKinnev fc Gammid (|2004 ): 
IKomissarovl (|2005|)). A generalization of the BZ process to 
bo osted non-spinning b lack holes has also been investigated 
by Neilsen et ajj ((20101), where the magnetic fields extract 
the translational kinetic energy from the black holes. In this 
case, there is also an extraction of rotational energy through 
the original BZ process if the boosted black holes are also 
spinning. During the coalescence of binary black hole sur- 
rounded by a magnetized circumbinary disk, this generalized 
BZ process will produce a dual jet structure during the in- 
spiral p hase which will result in to a single BZ jet after the 
merger (jPalenzuela et al.l ((20101)). 

The basic effects of the BZ mechanism can 
be explained by invokin g the membrane p ar adigm 
(see JThorne fc Macdonaldl (| 19821 ): JThorne et at] (|l986l ) 
for details), which endows to the black hole horizon some 
physical properties like a surface charge density and a resis- 
tivity. The problem is then reduced to a spherical conductor 
with a relative motion with respect to asymptotic magnetic 
field lines via rotation or translation. The magnetic field 
is produced by an external source and described by the 
force-free approximation. In spite of its simplicity and rel- 
ative success, this analogy does not yet explain the source 
of the energy, which can not be assigned to the hor i zon d ue 
to causality arguments ((Punslv fc Coronitil (| 19891 . (l990)). 
The membrane paradigm implies that the key ingredient 
of the mechanism is the black hole horizon, in contrast 
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with the arguments, pointing rather to the ergosphere, 
presented by Komissarov. Because the intrinsic marriage of 
the horizon and the ergosphere on black hole spacetimes, 
one could confuse the physical phenomena generated by 
each of them. It is therefore desirable to study the effect of 
each component separately. 

In this paper, we perform a systematic study of the iso- 
lated effect of the ergosphere in the EM fields, by consider- 
ing regular spacetimes produced by rapidly rotating neutron 
stars. By increasing the compactness of the star an ergo- 
sphere appears with a toroidal topology (see lAnsorg et al.l 
l|2002h for details). The compact object is immersed in a 
force-free environment produced by an externally sourced 
magnetic field. We will assume that the force-free EM fields 
are not coupled to the fluid, so their dynamics will be de- 
termined only by their evolution equations and by the prop- 
erties of the curved spacetime. Our aim is to analyze the 
precise role of the ergosphere on the activation of the BZ 
mechanism. 

The paper is organized as follows. A detailed descrip- 
tion of our model and a summary of the force-free evolution 
equations on a curved background is presented in Section [21 
Some known results for stationary and axisymmetric space- 
times are summarized in Section [S] The numerical setup 
and the initial data is described in Section |4l while that 
our numerical results are discussed in Section [5] Finally, we 
summarize our conclusions in Section [6] The robustness of 
our solutions against several sources of error is studied in 
the Appendix. 



is denoted by 7ij. Coordinates defined on adjacent hyper- 
surfaces can be related through the lapse function a, that 
measures the proper time elapsed between both hypersur- 
faces, and the shift vector /3\ that controls how the spatial 
coordinates propagate from one hypersurface to the next. 
An observer moving along the normal direction to the hy- 
persurfaces (Eulerian observer) will have a coordinate speed 
given by — /?\ and will measure a proper time dr = adt. In 
terms of these quantities, it is possible to bring the metric 
of the spacetime into the form 

ds — gab dx'^ dx 

= -a dt^ + 7,j (dx' + P' dt\ (dx^ -f P^ dt\ . (1) 

Here, and in what follows, Latin indices from the beginning 
of the alphabet(a, &, c, • • • ) denote four-dimensional space- 
time quantities, whereas Latin indices from the middle of 
the alphabet {i,j, k,- ■ ■) are spatial. It is also convenient to 
introduce the extrinsic curvature Kij, which is associated 
to the way in which the hypersurfaces are immersed in the 
spacetime {M,gab), in the form 

K^j = -:^(a.-£;3)7... (2) 

Za 

Notice that the Eulerian observer is defined independently 
of the space coordinates. It can be interpreted as being 
at rest in the hypersurface E^. In the context of spinning 
stars or black holes, this observer is also called the locally 
non-rotating observer or zero-angular-momentum observer 
(ZAMO). 



2 MODEL OF PASSIVE FORCE-FREE 
ENVIRONMENT 

We consider the evolution of a magnetized plasma with neg- 
ligible inertia on the spacetime produced by a rotating com- 
pact star which is assumed to be both stationary and ax- 
isymmetric. Our approach will involve the resolution of two 
different systems of equations. On one hand, the initial data 
is obtained by solving the Einstein equations coupled to the 
Hydrodynam ic equations. We wil l use an initial data solver 
developed bv lAnsorg et al.l (12003) in order to obtain the so- 
lution for both the fluid and the spacetime geometry. On 
the other hand, we will evolve the hyperbolic PDE system 
for the low-inertia magnetized plasma on this curved back- 
ground, which can be described by the force-free approxima- 
tion of the Maxwell equations. An important point of our 
model is that it neglects any coupling between the plasma 
and the fluid of the star. In this way, the electromagnetic 
fields will not interact directly with the fiuid, and its evo- 
lution will be determined solely by the force-free equations 
in a curved spacetime. In this section, we summarize the 
formulation used to describe these systems of equations. In 
particular, we review in detail the Eulerian description of 
electrodynamics in the force-free approximation. 



2.1 The 3+1 decomposition 

We consider a spacetime (M, gab) which is foliated by a fam- 
ily of spacelike hypersurfaces Et parametrized by time func- 
tion t. The induced metric on these spatial hypersurfaces 



2.2 3+1 decomposition of the Maxwell Equations 

The covariant Maxwell equations are given by 



VbF" 



4nr . 



ViT"'' ^0, 



(3) 



where l'' is the 4-current and F"'', * F'"'' are the MaxweU 
and the Faraday tensor, respectively. In order to provide an 
Eulerian description of the above equations, it is convenient 
to introduce the electric and magnetic fields measured by 
those observers, namely 



E" 



F Ub ., 



B = F Ub . 



(4) 



where n" is the unit vector normal to the hypersurface Et. 
Notice that, if the electric and magnetic susceptibilities of 
the medium vanish, as in vacuum or in a highly ionized 
plasma, the Faraday tensor becomes the dual of the Maxwell 
tensor. In a similar way, we define the charge density and 
current as 



ra 
1 Ua , 



J a I a rb 

—-i-b J 



(5) 



where -L(I= S''b + n"" Ub is the projection operator onto the 
hypersurface Et. Using the previous definitions, the Maxwell 
equations can be rewritten in the form 

{dt-Cp)E' = ^'''D,{aBu) + aKE'-A-KaJ\ (6) 

{dt-Cfi)B' = -€''''D,(aEk)^aKB\ (7) 

DiE' = 47rg, D^B' = . (8) 

Here Di =1-1 Va is the covariant derivative associated with 
the spatial metric 'ytj and e'-'*' is the Levi-Civita tensor. 
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It is useful to introduce, for later convenience, the vector 
potential Ua which can be decomposed into 



$ = -Ua n" 



Aa =-L' Ut 



(9) 



In terms of this vector potential, the Maxwell tensor can be 
written down as 

Fab = ~2\IyJAb]- (10) 

On the other hand, the electromagnetic energy-momentum 
tensor. 



Tab — - — 

47r 



Fa" Fbc — - gab F" Fed 



(11) 



(12) 



can be decomposed in the form 

Tab = £ ria nt + 2 n(a St) + Sab , 

where £, Sa and Sab correspond to the local electromag- 
netic energy density, the momentum density (Poynting vec- 
tor) and the spatial stress tensor as measured by the Eule- 
rian observer. Finally, the local conservation of the energy- 
momentum tensor (|lip is given by 



VbT"" = 



-F"' h . 



(13) 



The key point about this discussion is that it has been 
formulated in terms of physical quantities measured by 
the Eulerian observer (ZAMO). In order to close the sys- 
tem of the Maxwell equations, where a relation between 
the EM fields and the electric current is required, one 
can use quantities measured by the Eulerian observer in 
the sa me way as in the special re l ativistic electrodyn am- 
ics (see lMacdonald fc Thorn3 ljl982l ): lKomissarovl (|2004 l for 
details) . 



2.3 Force-free approximation 

The force-free approximation is valid in magnetized plas- 
mas when the electromagnetic energy density £ dominates 
over matter energy density. It happens, for instance, in the 
magnetospheres of neutron stars or black holes, where the 
density of the plasma is so extremely low that even mod- 
erate magnetic field stresses will dominate over the fluid 
pressure gradients. In this limit, the stress-energy tensor of 
the plasma therefore satisfies 



Tab rpab . rjiab 
— -'fluid + J- EM 



rpab 

Jem ■ 



(14) 

The local conservation of this stress-energy 
tensor implies that t he Lorentz forc e van - 
ishes, F°-^ Jb ~ llGoldreich fc JuhanI (|l969l ): 
IR. D. Blandford and R. L. Znaiek.1 1 19771 1). This expression 
can be written, in terms of 3-1-1 quantities, as 



E' Ji^Q, 



qE' + e''''JjB,,=0. 



(15) 



Taking the scalar and the vector product between the mag- 
netic fleld B'' and the spatial projection of the Lorentz 
force (|15|l . we obtain 



E'Bi=0, 



r 



5^(^ii+Jl 



(16) 



where Jl and J]_ are the component of the current paral- 
lel and perpendicular to the magnetic field _B\ respectively. 
These are defined as 



Jl = f Bi B' 



J'l^- 



ijk 



Ej -Bfc . 



(17) 



The first relation in (|16p implies that the electric and mag- 
netic field must be perpendicular. The second one defines the 
current up to the parallel component Jy . Using the Ma:xwell 
equations, one can compute {dt — Cp){E^ Bi) = 0, which has 
to vanish due to (|15p . This condition imposes a constraint for 
J|l , which can be substituted i nto Eq. (1 1611 to c omplete the 
specification of the current (see lGruzinovl (|2007|) for details). 
We will use instead an alternative prescription to enforce 
the force-free conditions, which has been used successfully 
in pre v ious studies of force- free magnetospheres (.Spitkovskvl 
(|2006l ): |Palenzuelaet all (|2010l ')'). 



3 STATIONARY AND AXISYMMETRIC 
SPACETIMES 

In the previous section, we have summarized the Maxwell 
equations and the force-free approximation on a generic 
spacetime. Nevertheless, since we are interested on station- 
ary and axisymmetric spacetimes, one can consider a set of 
coordinates adapted to these symmetries. In these coordi- 
nates, the metr ic of t h e spa c etime can be broug ht into the 
standard form (| Lewis! l|l932l 'l: IPapapetrorn (| 19661 )) 



ds 



■ dt + g4,4> {d(t> — Lodt) + Qr 



■ dr^ + gee dO^ , 



(18) 



where the metric coefficients {a, u, g^r, i7ee,5<p0} depend 
only on r and Q. Note that this metric describes usual as- 
trophysical objects such as neutron star or black hole space- 
times. In p articular, the Kerr metric ca n be written in the 
above form (jBergamini fc Viaggiul (|200J)). As we have men- 
tioned before, the shift vector is related to the relative veloc- 
ity between the Eulerian observer and the stationary spatial 
coordinates. One can then interpret lo as the drag velocity 
of this observer with respect to the hypersurface Et. 

Using the decomposition of the Maxwell tensor in terms 
of the vector potential ((9]), the condition of axisymmetry and 
stationarity implies that the electric field is purely poloidal, 
E^ = 0. According to H16p . it follows that _B' is perpendicu- 
lar to the poloidal components of the magnetic field, so that 
one can rewrite E^ in the form 



E-" = e'.fc B' U" , 



(19) 



where (7* is an axial vector given by (|Komissarovl (1200J)) 



a 



(20) 



and x" = dcp is the axial killing vector of the spacetime. No- 
tice that, according to (|18p . the shift vector is /3'° = —ujx''- 
Therefore, one can interpret the velocity U^ as the ve- 
locity of the magnetic field relative to the Eulerian ob- 
server and fi as the angular velocity of the magnetic lines, 
which can be written in terms of t he M axwell tensor as 
IJR. P. Blandford and R. L. Znaiekl (|l977l l) 



D. 



Ftr 

Ft4, 



Fte 
Fe^ 



(21) 



It is also useful to calculate the scalar B^ — E^ which, 
using the electric field defined by (|19|l . takes the form 

j2 e;2^^2 ^ B'^ a'^ - ggy^ Bp {D, ~ bj)'^ , 



{B^ 



(22) 



where Bp — Br B^ ~\- Be -B* is the magnitude of the poloidal 
component of the magnetic field. This relation implies a 
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change of sign of this invariant in highly compact rotating 
spacetimes with large g^^jo? and a;. Notice that, in elec- 
trovacuum scenarios, the Maxwell equations imply that -B** 
vanishes. In this case, one can also assume that the magnetic 
field is generated by distant plasma of large inertia, which 
means that the resulting magnetosphere will reach a steady 
state when fi = 0. This implies that the invariant (|22p be- 
comes 



{B' ~E'^)o? 



Bl{c?^f) 



(23) 

Inside the ergosphere cP' — ff' < 0. Therefore, the change of 
the sign of the invariant is related, at least in electrovacuum, 
with the presence of an ergosphere. 

Finally, it is possible to define conserved quantities 
associated to the killing vectors of the spacetime ^^ = 
dt and x" = (^<t>- On o^^ hand, the red shifted en- 
ergy density, c orresponding to the killing vector ^°, i s 
defined as fsee IR. D. Blandford and R. L. Znaiek.1 l|l977t ): 
iMacdonald fc Thornd (|l982l )) 



£-i~T°' ^.aTih = a £^ + o; S" Xi 
with fiux of energy given by 



SI 



-r'^i'bCc 



aS^ + uj S^^ Xj 



(24) 



(25) 



On the other hand, the angular momentum density, associ- 
ated to the killing vector x*, is defined 



^x = -T" Xa nb = S' Xi: 



(26) 



smoothly. This is the case for the Kerr-Schild coordinates, 
where the energy flux density S^ can be written as 

fir 

o 

(33) 



S. 






2Mr 
AB'B'^sin^ 



^] sin 



47r 
At the horizon, where r — th and A = 0, it becomes 



s; 



2tt 



(B' 



f{nH-0.) sin^e. 



(34) 



where th = M + -/M^ 



is the radius of the horizon 



and Qh = a/{2M th) can be interpreted as its rotation 
frequency, which is just the rotation velocity of a Eulerian 
observer at the apparent horizon. This result implies that, if 
< n < Q,H and B^ ^ 0, then there is an outward directed 
energy flux at the horizon. Therefore, rotational energy is 
being extracted from the black hole due to the magnetic 
field lines. The use of Kerr-Schild coordinates allow for direct 
computations of the flux at the horizon without any specia l 
treatmen t as in J R. D. Bland ford and R. L. Znaiek.l l| 19771 ): 
Macdona ld fc Thor nc (1982). However, one message from 
this simple calculation is that energy comes out of the event 
horizon, which is forbidden at the classical level since the 
horizon is a null surface. The problem lies in the fact that 
the energy flux defined on other surfaces is not obviously 
positive definite. 



with a fiux of angular momentum given by 

6c I z ^ ^ij 



SZ rT-ibC I Z 

X = r -L (, Xc 



S'' Xj ■ (27) 

Since E^ vanishes, the poloidal fiux vector S\ satisfies 

S^ = -^(B'xOS^ (28) 

Using this condition and (|16p . it is straightforward to show 
that 



J c — i i O Y . 



(29) 



Therefore, both the flux of red shifted energy and the flux of 
angular momentum are transported along the poloidal field 
lines. 

According to Eq. (|25|l , the EM radiated energy crossing 
a spherical surface at a given radius is 



dtS = 2Ti 



-gSldO. 



(30) 



Note that, on a regular spacetime in Lewis-Papapetrou co- 
ordinates (|18p . the radiated energy flux density S'^ is given 

by 



Sl^-^B'-B^a^g,^. 

Ill 



(31) 



Moreover, on the case of a Kerr spacetim e in Lewis- 
Papapetrou coordinates (|Bergamini fc Viaggiul (1200J)), the 
above expression becomes 



S. 



'2^ 



B'B'^A 



(32) 



where A = r^ + a? — 2M r. Since the Kerr spacetime in 
these coordinates is singular at the horizon, it is convenient 
to transform to other coordinates that penetrate the horizon 



4 NUMERICAL SET UP 

4.1 Diagnostic Quantities 

To extract physical information, we monitor the rotation 
frequency of the magnetic field lines (|2H) which is constant 
along magneti c fields lines on axisymmetric and station- 
ary solutions (|R. D. Blandford and R. L. Znaiek.l l|l977h '). 
and the Newman- Penrose electromagnetic scalars {$0,^2}, 
which are computed by contra cting the Maxwell tensor with 
a suitable null tetrad (see e. 0. iTeukolskvl l| 19731 )). 



$0 



Fab 7 



$2 






(35) 



The total energy flux (luminosity) of electromagnetic 
waves, which accounts for the energy carried off by outgoing 
waves to infinity, is 



Li 



lim 



$2 — $2 



I'l'o 



<\' 



dO. , (36) 



where $2 a^nd <l>(f are the background scalars produced by 
the steady part of the solution, which vanish only at far 
distances from the electromagnetic sources. However, since 
we are considering for simplicity that the magnetic field is 
produced by a very distant external source, there will be a 
non-zero contribution to these background scalars induced 
by the asymptotically uniform magnetic field configuration. 
An isolated system with no incoming radiation satisfies $0 = 
^o . Moreover, far from the star is valid the assumption that 
the background is approximately the same for the incoming 



and outgoing waves, so that <I>2 



$0 • Combining these 



relations with the general form given by (|36|l . it is obtained 
the simplified formula 



L^„ = lim 



r^\<^2-'^ofdn, 



(37) 
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which has been used previously in several works, re- 
producinR succes s fully the e x pected analyt i cal re lations 
l|Palenzuela et al.l (|2009l . I2OI0I ): iNeilsen et all (|2010l )'l. No- 
tice that these expressions are equivalent to the radiated 
energy (|30p evaluated at spatial infinity. 



4.2 Numerical methods 

We will use a finite difference scheme on a regular Carte- 
sian grid to solve numerically the hyperbolic PDE system. 
To ensure sufiicient resolution in an efficient manner we em- 
ploy Adaptive Mesh Refinement (AMR) via the had com- 
putational infrastruc ture, that pr o vides di stributed , Berger- 
Oliger style AMR (|HAD-Teaml (|2002l ): iLieblingj (|2002l )) 
with full sub-cycling in time, together with an improved 
tre atment of art i ficial boundaries as it has been presented 
bv iLehner et al.l l|200q l. For these simulations, the refine- 
ment regions are fixed initially and not changed during the 
evolution (i.e. Fixed Mesh Refinement). 

The spatial discretization is performed by using a fourth 
order accurate scheme satisfying the Summation By Parts 
rule. The time evolution is performed through the method of 
lines using a third order accurate Runge-Kutta integration 
scheme with a Courant parameter of A = 0.25 such that 
At — 0.25 Aa; holds in each refinement level. 

Our numerical domain consists of a cubical region de- 
fined by the intervals x^ £ [32 M, 32 M] with 61 points in 
the coarsest grid. We employ a fixed mesh refinement con- 
figuration with 6 levels of refinement, each one covering 
half of the domain of the parent coarser level. The coars- 
est resolution employed is Aa; = 1.07 while the finest one is 
Ax = 0.017. The radius of the different stars in these units 
are described in Table [T] We have adopted maximally dis- 
sipative boundary conditions in our simulations, by setting 
to zero the (time derivative) of the electrovacuum incoming 
modes (jPalenzuela et al.l HoTJ)). 



4.3 Initial data 

The initial data for the spacetime geometry and the fiuid 
variables produced by rotating stars is obtained by solv- 
ing the Einstein and the Hydrodynamic equations with the 
assumptions of stationarity and axisymmetry. The rotat- 
ing star sol utions have be e n coii structed using the code de- 
veloped by lAnsorg et al.l (|2002l ') based on a multi-domain 
spectral- method for representing the metric functions. The 
use of a spectral code has shown to be necessary to achieve 
high accuracy in the case of a stiff e quation of state (e.g. for 
constant total mass-energy density (JBonazzola fc Schneideij 
[l97J)). We consider equilibrium solutions for a rigidly rotat- 
ing star with an equation of state for homogeneous matter 
with constant total mass-energy density, jj, = const. For the 
calculation we use two different line elements to describe the 
exterior and the interior of the star. The Lewis-Papapetrou 
line element (llSp that covers the exterior has the form 



ds — 



+e 



'^ (ujdt — d(j)) 

{dp" + da 



e^" dr + We. 

2a /J 2 , ,y2 



(38) 



The advantage of this line element is that allows the metric 
potential i/ to remain real inside the ergosphere. For the 



interior of the star, in the comoving frame of coordinate, 
the metric can be expressed as 

,2 2(7 ,,2 , -2(7 r 2fe / , 2 , ,,2\ 

ds — ~e dt + e e (dp + d^ ) 

+ {W^ + ri) d<l>^] . (39) 

The potential U can be expressed in terms of the lapse func- 
tion a, while rj is the so-called gra vitomagnetic potent ial as- 
sociated with the shift vector (see iMeinel et al.l (|2008l ) for a 
detailed description). Given the particular equation of state 
and using the conservation of the energy-momentum tensor 
for the fluid, we obtain inside the star 



u 
e exp 



dp 



P+P. 



oVo 



const . 



(40) 



Isobaric surfaces inside the star correspond to a constant 
value of Vb. At the surface, where pressure goes to zero, it 
is possible to compute the redshift of a photon emitted with 
zero angular momentum via 



-Vo 



1, 



(41) 



By changing the parameter Vb, the solution becomes more 
compact and may contain an ergosphere. We have con- 
structed several rotating stars, with different value of Vb 
and rotation frequency Q. We kept the ratio between po- 
lar and equatorial radius constant. For all the models the 
value of the dimensionless spin parameter is roughly con- 
stant, a = J/M^ ~ 0.9. The mass, radius and other param- 
eters of the solutions are displayed in Table [ij where all the 
solutions listed in the table contain an ergosphere. Our most 
compact star is close to the limit of maximum compactness 
M/R < 4/9 « 0.44 for this family of solutions. 

Stationary and asymptotically ffat configurations with 
an ergosphere but without an horizon have been proved to be 
unstable or marginally unstable under to scalar and electro- 
magnetic perturbations iJFriedmanl (|l978l )). For slowly rotat- 
ing relativistic stars, the time scale of the i nstability is shown 
to be l onger than the Hubble time (see e.g. Comins fc Schutj 
l| 19781) 1. It has been shown bv lCardoso et all (120071 ) that, for 
the extreme case of compactness M/R > 0.5 and angular 
momentum J > 0.4M'^, the instability timescales reaches 
0.1 seconds for an object with mass of IM. In our simula- 
tions, however, both the spacetime and the fiuid are station- 
ary and therefore this instability cannot be active. 

The initial data for the black hole is analytical. In the 
case of spin ning black holes, w e will use the Kerr-Schild 
coordinates (|Keril(|l963l ). see e.o. lD. Kramer fc Herltl (|l980l ) 
for more details), 



Qab ~ Tjab + 2 H la lb , 



(42) 



where rjab is the Minkowsky metric and the scalar function 
H and the null vector l*" are defined by 

rM 



H 



lb 



r^ + a^zyr^ 


' 




f rx + ay 


ry — ax 


z 
r 



(43) 
(44) 



The compact object, either a neutron star or a black hole, 
is immersed in the external magnetic field produced by a 
distant current loop. This magnetic field is nearly constant 
initially near the compact object. In addition, it is chosen 
to be aligned with the spin of the compact object, which is 
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initially oriented along the z-axis. Therefore, the EM fields 
are initially set to B'' = Bo z and _B' = throughout all 
the domain. The field strength Bo is irrelevant, since we are 
assuming that the force- free fields behaves like test fields (i.e. 
they do not modify the curvature of the spacetime), and it 
has been set to Bo = 0.01. Since we are not considering 
any coupling between the fluid and the force-free EM fields, 
the dynamics of the latter will be only influenced by the 
regular spacetime both inside and outside the star. For all 
the effects, there will be no direct interaction between the 
EM fields and the fiuid. 



5 RESULTS 

In this section we will describe the dynamics of the force- free 
fields evolving in the stationary spacetimes produced by very 
compact rotating objects with dimensionless spin parameter 
a ~ 0.9. We will concentrate on the cases with the presence 
of an ergosphere in the spacetime (see the Appendix for a 
discussion on the cases without ergosphere). We will also 
analyze the EM power (if any) emitted by the BZ process 
in these spacetimes, as well as the features of the EM fields 
after they have relaxed to the stationary solution. 

All our simulations display an initial transient, in which 
the magnetic field is dragged and twisted around the spin- 
ning spacetime and induces a poloidal electric field. At late 
times the EM field relaxes to a stationary state, which 
is displayed in figure [T] for a representative case of regu- 
lar spacetime with an ergosphere. This case corresponds to 
Vo = —1.20 (see Table[l|. For comparison purposes, we have 
additionally included the results for the black hole case, a 
spacetime with an ergosphere but also with an horizon which 
hides a singularity. Further information on the structure of 
the solutions can be inferred from the currents and charge 
density of these two cases, as displayed in fig. (2] 

In the case of regular spacetimes with an ergosphere, 
the magnetic flux near the star is initially expelled, present- 
ing large damped oscillations that relaxes after few periods. 
During this relaxation, which seems to be more relevant as 
the compactness of the star increases, there is an impor- 
tant isotropic emission of energy. When the stationary state 
is reached, all the magnetic flelds from the region occupied 
by the star are twisted in the same direction than its an- 
gular momentum (in the 2: > domain). The currents, in 
this case, are composed of an outflow external cylinder and 
an inflow inner one. There is also a current sheet where 
B^ < E^ in the intersection of the ergosphere with the equa- 
torial plaJiej_sirnilai_to_yieone^ in the black hole 
case IJKomissarovl (|2004l '): |PaIenzuela et all (|2010l ')). 

The black hole simulation, on the other hand, relaxes 
to the stationary state in a shorter timescale than the above 
case. The final state resembles the solution corresponding 
to the regular spacetime with an ergosphere, displaying an 
analogous structure of magnetic flelds, currents and charge 
densities. This clearly indicates that the BZ mechanism act- 
ing on the spacetimes with an ergosphere is basically the 
same than in the black hole case. 

The poloidal structure of the magnetic flelds is almost 
identical in all the simulations, showing that the magnetic 
flux threading the spacetime occupied by the compact object 
is basically the same. The luminosity, evaluated in a sphere 



located at _R ~ 10 r^ for the stars, and conveniently rescaled 
for the black hole, is displayed in the left panel of fig. [3] for 
all the simulations. The luminosity increases very fast as the 
compactness of the star increases, although it does not reach 
the high values of the black hole. 

This smoothness is also found in the angular velocity 
of the magnetic field, displayed in the right panel of fig. [3l 
where Q has been normalized with respect to the central 
maximum value uuc for the stars, and with respect to flH for 
the black hole. As it was mentioned before, the angular ve- 
locity SI is confined to a small cylinder, showing that the jet 
is coUimated to the region occupied by the compact object. 
The fast growth of the maximum of this quantity as a func- 
tion of the compactness of the star can be fitted accurately 
in this regime to an exponential function, as it is shown in 
the left panel of figure [l] The luminosity for the different 
cases can also be represented as a function of the compact- 
ness, showing roughly also an exponential dependence in the 
right panel of figure |4l 

From our numerical results we have found the following 
scaling relations for the angular velocity i}, 



Q/oJc 



Ae 



\M/R 



(45) 



and for the ratio of poloidal and toroidal components of the 
magnetic field 



5<F 



-/OB'- 



(46) 



with / ~ 1/5 for the spacetime with an ergosphere. Notice 
that these estimates contain large sources of error, since they 
both neglects the details of the spacetime geometry and the 
azimuthal dependence of these quantities. Nevertheless, they 
can be used to study the behavior of the solution in different 
limits and to obtain the right order of magnitude of the 
luminosity. By using the line element of our iiutial data (|38|) . 
the energy fiux density (|3ip reduces to 



-—B^'B'^W^ 

2ti 



fJL{B^fW^ 



2 TT 



(47) 



which is a positive definite quantity, and where we have 
used subsequently the above approximations. Notice that 
the scaling is similar to the BZ power in Eq. (|34p . except 
by the exponential dependence on the compactness, and is 
consistent with the results displayed at the right panel of 

fig-H 



6 CONCLUDING REMARKS 

We have studied the evolution of EM fields on rotating and 
highly compact regular spacetimes with an ergosphere. Our 
results show that if a ergosphere is present the structure of 
EM fields and currents are similar to the black hole ones. 
This implies that the same mechanism operates in both 
spacetimes, independently of the presence/absence of an 
horizon. Notice that these results are in agreement with the 
fact that the BZ process is not an effect caused by the h ori- 
zon, as it was pointed out by iKomissarovl l|2004l . [2005|) in 
the context of black holes. 

In the case of a realistic rotating star, the fluid will be 
coupled with the EM fields and, therefore, they are forced 
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Figure 1. Rotating star and black hole. Magnetic field lines component in the plane a; = after the quasi-stationary state is reached, 
corresponding to the spacetiine with an ergosphere (left panel) and to the black hole (right panel). The vertical lines indicate the poloidal 
component, while that the blue-red colors indicates the strength of the component normal to the plane. The structure of these two 
components of the magnetic field for both cases, star and black hole, are quite similar each other. The surface of the star is plotted as a 
black ellipsoid, while that the ergosphere is plotted in red. 
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Figure 2. Rotating stars and black hole. Induced charge density (in red-blue colors) and the poloidal currents (in vectors) on the plane 
X = Q a,t the quasi-stationary state, corresponding to the same cases than figure [T] The surface of the star is plotted in black and the 
ergosphere in red lines. 




20 30 

time 



-BH 


1=0.9 


- Vo= 


-1.20 


- Vo= 


-1.1.5 


- Vo= 


-1.10 


Vo= 


-1.0.5 


- Vo= 


-1.00 



9 (radians) 



Figure 3. Left panel displays the EM luminosity obtained in the rotating spacetimes with a = J/M^ ?» 0.9. The luminosity increases 
monotonically with the compactness. Right panel displays the angular velocity of the magnetic field Q, computed at r Ri Sr^ and 
normalized with respect to its maximum value inside the star (see Table [TJ . 



© 2002 RAS, MNRAS 000, [lH9] 



8 Milton Ruiz, Carlos Palenzuela, Filippo Galeazzi, and Carles Bona 

Table 1. Afj, baryonic mass, Madai gravitational mass, J angular momentum, Q angular velocity of the 
fluid, re equatorial radius, re circumferential radius in Schwarzschild coordinates, Wc metric potential and 
Vb parameter. All the configurations contain an ergosphere.The ratio between polar and equatorial radius, 
Tp/re , is equal to 0.49 for all the configurations. All numerical values are in units with G = c = 1. 



Vo 
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Madm 
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re 


re 
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Figure 4. (Left panel) The maximum of the (normalized) angular velocity of the magnetic field U/lUc, as a function of the compactness. 
A fit with the numerical results shows an exponential dependence Ae^^^'^ with A Ri 3 X 10^^" and A fs 50. (Right panel) The 
EM luminosity calculated in the regular spacetimes, normalized with respect to the black hole one. Roughly, it seems to depend also 
exponentially on the compactness. 



to rotate following the fluid. So, the extraction process can 
be described by the simple model of the Faraday disk, which 
is even more efficient than the BZ process. In this case, the 
twist on the magnetic field will not be dominated by the 
spacetime but by the fluid rotation. In general, the strong 
dependence of Q on the compactness will imply a low en- 
ergy extraction on regular spacetimes even if matter is not 
present. 

An example of regular and rotating spacetime may be 
produced by an orbiting bina ry black hole system, as it was 
considered bv lLvutikovl (|201ll ). In this work it was suggested 
that the rotation of the spacetime inside the binary system, 
induced by the orbital motion of the black holes, would allow 
for two different channels of energy extraction; through the 
generalized BZ in boosted (and maybe rota ting) black holes, 
which will produce a dual jet structure (Palenzue la et al.l 
l|201Cl )). and through the BZ process on regular spacetimes 
inside the binary system. However, the lack of an ergosphere 
in the central region will probably prevent the activation of 
the BZ process (see the Appendix for a discussion on the case 
of spacetimes without ergosphere). Even if the BZ process 
takes place, and assuming that the luminosity still follows 
the exponential dependence shown earlier, the low compact- 
ness of the spacetime in this central region will induce a faint 
jet that will be overshined by the BZ process taping kinetic 
energy (either translational or rotational) directly from the 
black holes. 

Summarizing, we have found that the BZ process is 
also present on regular spacetimes with an ergosphere. Our 
conclusive results implies that we have to reconsider the 
membrane paradigm as a tool to explain the BZ mecha- 



nism which also seems to be able to extract energy from 
rotating regular spacetimes with ergospheres and boosted 
non-rotating black holes. 



ACKNOWLEDGMENTS 

We are in debt with L. Lehner and T. Garrett for stimu- 
lating discussions about the origin of the BZ process which 
initiated this project. We are also grateful to Marcus Ansorg 
and his collaborators for giving us the access to his spectral 
code and for the useful discussions on spectral method dur- 
ing his period at AEI. It is a pleasure to thank L. Lehner, 
S. Liebling, S. Komissarov, C. Thompson, as well as to our 
long term collaborators M. Anderson, E. Hirschmann and D. 
Neilsen for useful discussions and comments on this paper. 
We acknowledge support from Spanish Ministry of Science 
and Innovation under grants CSD2007-00042, CSD2009- 
00064 and FPA2010-16495, and Govern de les Illes Balears. 
Computations were performed in Scinet and Mare Nostrum 
(funded by BSC Grant No. FI-2011-3-0017). 



REFERENCES 

Ansorg M., Kleinwachter A., Meinel R., 2002, Astron. As- 
trophys., 381, L49 

Bergamini R., Viaggiu S., 2004, Class. Quant. Grav., 21, 
4567 

Bonazzola S., Schneider S., 1974, The Astrophysical Jour- 
nal, 191, 273 



© 2002 RAS, MNRAS 000, [l]l] 



The role of the ergosphere in the Blandford-Znajek process 9 



Cardoso V., Pani P., Cadoni M., Cavaglia M., 2007, Phys- 
ical Review D, 77, 14 

Comins N., Schutz B., 1978, Proceedings Of The Royal So- 
ciety Of London A Mathematical And Physical Sciences, 
364, 211 

D. Kramer, H. Stephani M. M., Herlt E., 1980, Exact Solu- 
tions of Einstein's Field Equations. Cambridge University 
Press, Cambridge 

Friedman J. L., 1978, Communications in Mathematical 
Physics, 63, 243 

Goldreich P., Juhan W. H., 1969, Astrophys. J., 157, 869 

Gruzinov A., 2007 

HAD-Team, 2002, The had code 

Kerr R. P., 1963, Phys. Rev. Lett., 11, 237 

Komissarov S. S., 2002, ArXiv Astrophysics e-prints 

Komissarov S. S., 2004, Mon. Not. Roy. Astron. Soc, 350, 
407 

Komissarov S. S., 2005, Mon. Not. Roy. Astron. Soc, 359, 
801 

Komissarov S. S., 2009, Journal of Korean Physical Society, 
54, 2503 

Lehner L., Liebling S. L., Reula O., 2006, Class. Quant. 
Grav., 23, S421 

Lewis T., 1932, Proc. Roy. Soc. London Ser. A, 136, 176 

Liebhng S. L., 2002, Phys. Rev., D66, 041703 

Lyutikov M., 2011, Phys. Rev., D83, 064001 

Macdonald D., Thorne K. S., 1982, Mon. Not. Roy. Astron. 
Soc, 198, 345 

McKinney J. C, Gammie C. F., 2004, The Astrophysical 
Journal, 611, 977 

Meinel R., Ansorg M., Kleinwachter A., Neugebauer G., 
Petroff D., 2008, Relativistic Figures of Equilibrium, 
Meinel, R., Ansorg, M., Kleinwachter, A., Neugebauer, 
G., & Petroff, D., ed. 

Neilsen D., Lehner L., Palenzuela C, Hirschmann E. W., 
Liebling S. L., Motl P. M., Garrett T., 2010, ArXiv e- 
prints 

Palenzuela C, Anderson M., Lehner L., Liebling S. L., 
Neilsen D., 2009, Phys. Rev. Lett., 103, 081101 

Palenzuela C, Bona C, Lehner L., Reula O., 2011, Classi- 
cal and Quantum Gravity, 28, 134007 

Palenzuela C, Garrett T., Lehner L., Liebling S. L., 2010, 
Phys. Rev., D82, 044045 

Palenzuela C, Lehner L., Liebling S. L., 2010, Science, 329, 
927 

Papapetrou A., 1966, Ann. Inst. H. Poincare A, 4, 83 

Punsly B., Coroniti F. V., 1989, Phys. Rev., D40, 3834 

Punsly B., Coroniti F. V., 1990, The Astrophysical Journal, 
350, 518 

R. D. Blandford and R. L. Znajek., 1977, Mon. Not. R. 
Astr. Soc, 179, 433 

Spitkovsky A., 2006, The Astrophysical Journal, 648, L51 

Tchekhovskoy A., Narayan R., McKinney J. C, 2010, The 
Astrophysical Journal, 711, 50 

Teukolsky S. A., 1973, Astrophys. J., 185, 635 

Thorne K. S., Macdonald D., 1982, Mon. Not. R. Astr. 
Soc, 198, 339 

Thorne K. S., Price R. H., Macdonald D. A., 1986, The 
Membrane Paradigm. Yale University Press, New Haven 



APPENDIX A: 
RESULTS 



ROBUSTNESS OF THE 



We have checked the robustness of our solutions against sev- 
eral sources of error. One of the possible problems may come 
from the way in which the analysis quantities are evaluated. 
In particular, the surface where the luminosity is computed 
may be located too close to the source, where the spacetime 
is still far from being flat, producing an error in that mea- 
sure. We have compared the luminosity computed in two sur- 
faces located and 5 re and 10 re for the case with Vo = —1.2, 
obtaining a difference smaller than 5%. Another potential 
problem may come from the influence of the boundary condi- 
tions, which may produce unphysical reflections which may 
affect, after a light crossing time, the dynamics of the sys- 
tem. In our simulations this is not a problem since the so- 
lution relaxes to the stationary state before a light crossing 
time, and it remains unaffected afterwards. 

Probably the most important source of inaccuracies 
comes from numerical discretization errors. We have com- 
pared three different spatial resolutions, corresponding to 
A'^ — {40, 60, 80} points in the coarsest grid. Our compar- 
isons are summarized in fig. lAll where we have restricted our 
analysis to a representative case corresponding to Vb = —1.2. 
Notice that, in this case, there is an ergosphere, and conse- 
quently, a current sheet on the equatorial plane which is 
difficult to represent on a discretized grid. Nevertheless, the 
luminosity display the expected fourth order convergence to 
a well defined solution. 

We have also tried to study the relaxed solutions of 
spacetimes without an ergosphere. Although the luminos- 
ity reaches a quasi- stationary value, it changes dramatically 
with resolution and it does not seem to converge to a unique 
solution. This lack of convergence is shown in the left panel 
of fig. IA2I for a representative case without ergosphere cor- 
responding to Vb = —0.97. In these cases there is a violation 
of the force-free condition, as it can be seen in the right 
panel of figure IA2I The poloidal currents inside the star are 
not parallel to the poloidal magnetic field, a consequence 
of a non-vanishing toroidal electric field. This component 
has to vanish in stationary and axisymmetric spacetimes. 
Our guess is that there is no force-free solution with these 
boundary conditions unless a current sheet appears, and so 
the numerical evolution relaxes only to an approximated so- 
lution which depends strongly on the resolution. More work 
is needed to elucidate this issue, maybe considering the full 
MHD problem instead of the force- free limit. 

This paper has been typeset from a Tff)i/ M^^X file prepared 
by the author. 
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Figure Al. The EM luminosity for the cases Vb = —1.2 (with ergosphere) for three different spatial numerical resolution with 
N = {40, 60, 80} points in the coarsest grid. The luminosity converges with the expected fourth order convergence 
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Figure A2. (Left panel) EM luminosity for the same three different resolutions as in fig. lAll In this case the solutions do not display 
any convergence. (Right panel) Poloidal magnetic field lines (in blue) and poloidal currents vectors (in black) in the plane x = after the 
quasi-stationary state is reached, corresponding to the spacetime without an ergosphere. The blue-red colors indicates the charge density, 
while that the black ellipsoid represents the surface of the star. Although the poloidal component of these fields has to be parallel in a 
force-free stationary and axisymmetric solution, it is clearly not satisfied in the interior of the star. 
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